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Abstract 

We suggest a method to create turbulence in a Bose-Einstein condensate. The method consists 
in, firstly, creating an ordered vortex array, and, secondly, imprinting a phase difference in different 
regions of the condensate. By solving numerically the two-dimensional Gross-Pitaevskii equation 
we show that the motion of the resulting positive and negative vortices is disordered. 

PACS numbers: 
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I. INTRODUCTION 



The study of the turbulence in superfTuids - atomic Bose-Einstein condensates and liquid 



helium II |l( - help us to understand better issues of classical Euler fluid dynamics |2|. In 
a superffuid the rotational motion can only take the form of vortex filaments with fixed 
quantized circulation. The length of these filaments is macroscopic (they can be as long 
as the experimental cell) but the thickness of the vortex core is microscopic (about 10 -8 
cm in liquid helium II). In 3D superfluid turbulence is manifested as a disordered tangle of 
vortex filaments. In 2D disordered motion of the vortices is a sign of turbulence. Indeed, 
many remarkable similarities between classical turbulence and superfluid turbulence have 
been noticed. For example the classical k~ 5 ^ 3 Kolmogorov energy spectrum (where k is the 
wavenumber) was observed in helium II when agitated by rotating propellers [3]. The same 



spectrum was a 
collaborators 



3 



aarent in the numerical simulations of Nore et al [4] and of Tsubota and 
Another example is the classical t~ 3 ^ 2 temporal decay of turbulence 



(where t is time) which was measured in helium II behind a towed grid [3] 

In atomic Bose-Einstein condensates, it has been suggested that turbulence can be ere- 

- * — , „- .... - ™ a. - — — . «. ff 

the interactions of vortices with rarefaction solitary waves was investigated by Berloff [9j. 



The decay and forced 2D turbulent flow [10( was also described with numerical simulations. 
Large scale turbulence of quantized vortices was studied in superfluid 3 He-B and 4 He. The 
disadvantage of atomic Bose-Einstein condensate is that with existing techniques, conden- 
sates are rather small and contain at the most some hundreds vortices. On the other hand 
good visualization of individual vortices is possible. An interesting problem to study is the 
;ransformation of kinetic energy into acustic energy [4] arising from vortex reconnections 



ll| and vortex acceleration 12J. 



The aim of this paper is to show that the phase imprinting method Q] can be used 
to study the formation of rotating turbulence accompanied by acoustic emissions in a simple 
2-dimensional condensate. 
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II. MODEL AND NUMERICAL TECHNIQUE 

To explore the density and phase profile of a rotating condensate, we solve numerically 
the following two-dimensional Gross-Pitaevskii equation which governs the time evolution 
of the complex order parameter ip(r,t): 



h 2 

■—V 2 + V tr + g\m\ 2 -^-VLL z 
2m 



j *, (1) 
where the coupling constant is g — 4-Kh 2 a/m, m is the atomic mass and a is the scattering 



length. The chemical potential, /x, which can be initially estimated in dimensional form as: 



fjL = 2HJ-, (2) 
V m 

from the Thomas- Fermi solution [l8j] of Eq. [TJ is found by the single-step integration: 



In 



(\m\ 2 ) 



after the initial Thomas-Fermi solution has relaxed to a time - independent solution in the 
harmonic trap, where (...) denotes spatial average. The trapping potential is 

Vtr(x,y) = -ttujI {x 2 + y 2 ) , (4) 

where u± is the radial trap frequency. The quantity Q is the angular frequency of rotation 
about the z-axis and L z = ih(xd y — yd x ) is the angular momentum operator. 

In order to perform our analysis, we introduce the density and the velocity: p(r, t) = 
\ip(r, t)\ 2 and v = h/m'VS(r,t) and decompose the total energy, E tot , in the following way: 

Etot = E kin + E int + E q + E tr . (5) 

where the kinetic energy, Euni the internal energy, E int , the quantum energy, E q , and the 
trap energy, E tr , are given respectively by 

/h 2 / \ 2 

— ( V / P(x,t)v(x,t)j d 2 r, (6) 
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E mt (t) = J g(p(x,t)) 2 d\ (7) 

Eqit) = /l^( Vv/ ^) 2d2r ' (8) 

Etr{t) = J p(x,t)V tr (*)d 2 r. (9) 

Furthermore, we decompose the kinetic energy, E kin , into a part due to the sound field, 
Esound, and a part due to vortices, E vortex 



Ekin E soun( i -\- E vor f ex (10) 

At a given time t, the vortex energy, E vortex , is obtained by relaxing the Gross-Pitaevskii 
equation in imaginary time, which yields the lowest energy state for a given vortex configu- 
191 ]. The sound energy is then recovered from E sound = E kin — E vortex . 



ration 



It is convenient to rewrite the Gross-Pitaevskii equation in dimensionless form in terms 
of the harmonic oscillator energy hu>±, the harmonic oscillator length *Jh] (2mu±) and the 
harmonic oscillator time uT 1 . We obtain: 



.dip 



-^V 2 + V tr + C\ifj\ 2 -p-QL z 



^ (11) 



where 



v,r(x,y) = \(x* + y*), (12) 
C = ^, (13) 

N is the number of atoms and L the extension of the condensate in the z - direction. Unless 
stated otherwise, we set C = 1400 in our calculations. In dimensionless units the initial 
estimate for \i is ji = 1/2^J AC . The dimensionless form of the energy contributors are: 

E k in{t) = y^( v / p(x,t)v(x,t)) 2 rf 2 r, (14) 
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E int {t) = J C(p(x,t)) 2 d 2 r, (15) 

E q (t)= j\{v^^)) 2 d\ (16) 

and 

Etr(t) = J p(x,t)V tr (x)d 2 r. (17) 

The calculation is performed in a square box of size D. We choose D so that it is larger than 
the trapped condensate, and impose boundary conditions ip — at x — ±D/2, y = ±D/2. 
Typically, D = 10. The time stepping is a semi- implicit the Crank-Nicholson algorithm. 



III. RESULTS 



Following Leadbeater et al., 11], starting from a non-rotating condensate at time t = 
0, we set Q = 0.85 and create a stable lattice of 20 vortices in a rotating frame. This 
state mimics the solid body rotation of a rotating classical Euler fluid. The relaxation of 
the vortex system into a stable lattice is helped by the addition of a tiny dissipative term, 
— T^ Jf^ with 7 1 at the left hand side of Eq. [1] and Eq. [TTJ, as described by Tsubota et 



al, 151 ] ; this term models the interaction of the condensate with the surrounding thermal 



cloud at finite temperatures as described by Madarassy and Barenghi, [16j. Density and 
phase profiles of the condensate are shown in Fig. [TJ At t = 200, we suddenly perform the 
following phase imprinting: 

4* ~* 4> f or x<0,y<0 and x>0,y>0, (18) 

■0 — > e™^ for x < 0, y > and x > 0, y < 0. (19) 

For a large system with many vortices (as our system), the phase imprinting denotes that 
about half of quantized vortices change their sense of rotation. Fig. [2] shows the condensate 
just after the phase imprinting. Note that, the original condensate is divided into four equal 
parts, two parts rotating clockwise and two anticlockwise. This implies that the kinetic 
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energy of the system is suddenly lowered. Hereafter, we shall consider the contributions to 
the energy for t > 200, after the phase imprinting. 

After the phase imprinting, positive and negative vortices interact and move due the 
velocity field which they induce on each other. Fig. [3] shows phase profiles at t = 200.65 and 
t = 200.9. Fig. H] shows density and phase profiles at t = 201.95. Fig. Hk. reveals that large 
density oscillations (sound waves) are present in the condensate. 

Due to the confined and finite size of the condensate, the sound waves are reflected from 
the edge of the condensate and reinteract with the vortex pairs. We can see that when 
Esound is approaching its maximum, E vortex is approaching its minimum, see Fig. O The 
interaction of positive and negative vortices leads to a disordered motion of the vortices. 
The trajectories of the vortices can be obtained by tracking the phase profiles in time. For 
clarity, we show vortex trajectories of two vortices at the time in Fig. El Fig. [7] and Fig. [HJ 
Fig. [9] shows all these trajectories together. 



IV. DISCUSSION 



Our calculation shows that the method which we suggest can be used to create a system of 
positive and negative vortices in a disk-shaped condensate and that the vortices will move 
in irregular way under each other's influence. The method can thus be used to create a 
two-dimensional turbulent system, similar to the classical Onsager vortex gas 171 ] . Clearly, 
a disk-shaped condensate created in the laboratory is a three-dimensional system, but, 
provided that, the thickness of the disk is only few times the coherence length, the resulting 
turbulence will be sufficiently two-dimensional, because Kelvin waves will not be able to 
propagate along vortices in the z-direction. 

It is useful to distinguish two possible evolutions after the phase imprinting: 

(i): Q and 7 are kept the same. After a transient turbulent state, the system will relax 
to the original vortex lattice as in Fig. [1] (essentialy repeating the original calculation but 
starting from a different initial condition). Negative vortices will move to the edge of the 
condensate and vanish, as shown in Fig. (TDJ 
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(ii): Q is set to zero together with the dissipation, 7. If 7 ^ 0, the system ought to 
eventually relax to the ground state with no vortices. If 7 = 0, the negative vortices remain 
in the condensate, as in Fig. [TTj 



It is instructive to compare the number of total, negative and positive vortices in these 
two cases (Q and 7 kept the same and Q and 7 set to zero) for the same interval of time, 
between t = 199 and t = 205. We conclude, that in the case of (i), the negative vortices 
disappear quickly, but in the case of (ii), for that short time, the negative vortices stay in the 
condensate and the transient turbulent state lasts much longer. Over a much longer time 
scale, we expect that the condensate will become vortex-free, due to collisions and mutual 
annihilation of vortices with anti vortices. 

There are many possible variations of the method. For example, we can obtain positive 
and negative vortices by setting: ip — > tp for y > and ip — > ip* for y < 0, (where the star 
denotes the complex conjugate). If we set: ip — > ip for y > and ip — > e^ip for y < a dark 



soliton-like disturbance 



201 ] is created along the x-axis. In a two-dimensional system this 



structure is unstable to the snake instability [21] and vortex - anti vortex pairs (or vortex 



rings in three dimensions 22| are created. 
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20 40 60 80 20 40 60 80 

FIG. 1: Stable condensate with N = 20 vortices for f2 = 0.85 at t = 199.8 before the phase 
imprinting. The vortices are visible as holes in the density profile (left) or as discontinuities from 
to 2tt of the phase profile (right). 




FIG. 2: Density (left) and phase (right) profiles at t = 200.4. 
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FIG. 3: Phase profile at t = 200.65 (left) and 200.9 (right). 




FIG. 4: Density profile (left) and phase profile (right) at t = 201.95. 



11 




hi 1.5 - 

1 - 

0.5 -; ; 

----- t ---- — 1 

200 250 300 350 400 
t 

FIG. 5: Kinetic energy (solid curve, top), vortex energy (dashed curve, middle) and sound energy 
(dotted curve, bottom). 
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FIG. 6: Left: Trajectories of positive vortex a (solid curve) initially located at x = -1.25, y = 2.5 
and of negative vortex b (dashed curve) initially located at x = -1.12, y = 1.31. Right: Trajectories 
of positive vortex c (solid curve) initially located at x = -2.31, y = and of positive vortex d 
(dashed curve) initially located at x = 0, y = -1.27, from t = 200.5 to t = 204.7. 
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FIG. 7: Left: Trajectories of positive vortex e (solid curve) initially located at x = 1.38 , y = 0.12 
and of positive vortex / (dashed curve) initially located at x = 3.88, y = 2.45. Right: Trajectories 
of positive vortex g (solid curve) initially located at x = -2.35, y = -1.32 and of negative vortex h 
(dashed curve) initially located at x = -0.43, y = -1.17, from t = 200.5 to t = 204.7. 




i 1 1 1 1 r 




8 
6 
4 
2 

>- - 
-2 
-4 
-6 
-8 

-8 -6 -4 -2 2 4 6 8 

x 

FIG. 8: Left: Trajectories of positive vortex i (solid curve) initially located at x = 0.43 , y = -0.38 
and of negative vortex j (dashed curve) initially located at x = 0.42, y = -2.42. Right: Trajectories 
of positive vortex k (solid curve) initially located at x = 1.18, y = -4.85 and of negative vortex / 
(dashed curve) initially located at x = 1.35, y = -3.84, from t = 200.5 to t = 204.7. 
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FIG. 9: Vortex trajectories of vortices a, b, c, d, e, /, g, h, i, j, k, and I, for times: 200.5 < t < 
204.7. 
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FIG. 10: Number of total vortices, N (circles), positive vortices, N + (triangles) and negative 
vortices, N~ (filled circles) as a function of time if, after the phase imprinting at t = 200, O is not 
changed and remains = 0.85. 
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FIG. 11: The number of total vortices, N (circles), positive vortices, N + (triangles) and negative 
vortices, N~ (filled circles) as a function of time if, after the phase imprinting at t = 200, Q is set 
to Q = 0. 
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